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RESUME 

L'6tude  du  changement  du  coefficient  d' extinction  des  aerosols 
infrarouges  en  fonction  de  l'humidite  relative  est  l'objet  de  ce  rapport. 
On  se  base  sur  un  modele  a deux  types  de  distribution  des  particules: 
maritime  et  continental.  Le  facteur  du  coefficient  d' extinction  des 
aerorols  s'accroit  de  plus  de  10  dans  les  bandes  spectrales  de  3-5  et 
8-15  ym  lorsque  l'humidite  relative  passe  de  0.5  a 0.99.  Aux  humidites 
relatives  tres  elevees,  une  region  d' extinction  plus  faible  est  observee 
autour  de  10-11  ym.  (NC) 


ABSTRACT 

Y 

The  variation  of  the  infrared  aerosol  extinction  coefficient 
as  a function  of  relative  humidity  is  studied.  This  is  based  on  a two- 
component  particle  distribution  model.  The  aerosol  extinction  coefficient 
increases  by  a factor  of  more  than  10  throughout  the  3-5  and  8-15  t-.  >'•  'tii-'-/ 

spectral  regions  as  the  relative  humidity  increases  from  0.5  to  0.99. 

For  high  relative  humidities,  a region  of  lower  extinction  around 
10-11  w*r  is  noticed.  (U) 
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1.0  INTRODUCTION 

The  attenuation  of  an  electromagnetic  wave  propagating  in 
the  atmosphere  arises  from  molecular  and  aerosol  extinctions.  In  the 
3-5  pm  and  8-15  pm  'window'  regions  where  no  major  molecular  resonance 
absorption  is  present,  aerosol  extinction  is  of  primary  importance. 

With  a change  of  relative  humidity  in  the  atmosphere,  conden- 
sation or  evaporation  of  water  takes  place  on  the  aerosol  particles 
and  these,  at  the  same  time,  change  the  particle  optical  parameters 
which  determine  the  aerosol  extinction.  When  the  relative  humidity 
increases,  the  size  of  an  aerosol  particle  increases  through  the 
accretion  of  water.  A particle  can  triple  or  quadruple  its  dry-state 
size  at  relative  humidities  above  95%.  Together  with  the  growth  in 
size,  the  complex  refractive  index  of  the  particle  varies.  Using  a 
two-component  distribution  model,  Barnhardt  and  Streete  [1]  have 
calculated  the  aerosol  scattering  coefficient  in  the  infrared  as  a 
function  of  relative  humidity.  This  work  has  been  extended  by 
Hodges  [2]  to  include  the  absorption  of  the  aerosols.  However,  a 
basic  assumption  made  in  both  works  is  that  aerosol  particles  grow  at 
a rate  dependent  only  on  the  relative  humidity  and  not  on  the  chemical 
composition  of  the  particles  in  the  original  dry  state.  A recent 
comprehensive  study  has  been  carried  out  by  HSnel  [3]  providing 
empirical  data  on  the  change  of  aerosol  optical  parameters  as  functions 
of  relative  humidity  through  the  concept  of  the  mass  increase  coef- 
ficient. By  means  of  these  empirical  results,  a more  realistic  assess- 
ment of  the  effects  of  relative  humidity  on  aerosol  extinction  can 
now  be  made. 

In  practice,  the  aerosol  extinction  coefficients  over  the 
3-5  pm  and  8-15  pm  spectral  regions  are  of  particular  importance  to 
the  operation  of  broad -band  infrared  systems.  In  this  report,  we 
present  the  aerosol  extinction  coefficients  as  functions  of  relative 
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humidity  in  the  above  two  window  regions  obtained,  using  the  two-component 
particle  number  distribution  model  of  [1]  and  the  empirical  data  of  the 
mass  increase  coefficients  of  [3].  This  enables  us  to  take  into  account 
the  different  behavior  of  the  maritime  and  continental  aerosols  with 
the  change  of  the  relative  humidity. 

The  work  was  performed  at  DREV  between  January  and  November 
1977  under  PCN  33A11  "Studies  of  Aerosols". 


2.0  AEROSOL  EXTINCTION  AND  RELATIVE  HUMIDITY 


Let  n(f,  r)  be  the  number  of  aerosol  particles  per  unit  volume 
with  radius  between  r and  r + dr  at  a given  relative  humidity  f.  The 
corresponding  aerosol  extinction  coefficient  B(f,  X)  for  electromagnetic 
radiation  of  wavelength  X is  given  by 


3(f,  X)  = 7i 


"(f.  r)  Qext 


n(f.  X) 


] 


r2dr 


(1) 


where  a = 27rr/X  is  the  size  parameter;  n(f,  X),  the  complex  refractive 
index  and  Qext»  the  efficiency  factor  of  extinction. 


We  wish  to  consider  the  problem  by  starting  from  the  situation 
with  f = 0* . The  dry-state  condensation  nuclei  distribution  is 
n(f  = 0,  r).  To  find  the  extinction  coefficient  at  any  other  relative 
humidity  f > 0,  it  is  necessary  to  relate  n(f,  r)  to  n(0,  r)  and 
n(f,  X)  to  n(0,  X).  To  concentrate  our  attention  on  the  effects  of 
the  change  of  relative  humidity,  we  assume  that  as  f increases 
n(f,  r)  changes  only  by  virtue  of  the  condensation  of  water  vapor. 

This  is  to  say  that  all  the  particles  retain  their  individuality 
without  coagulation  and  furthermore,  no  new  particles  are  produced 
in  the  meantime. 

* We  have  formulated  our  problem  in  such  a manner  as  to  exclude  the 
consideration  of  the  hysteresis  effect.  8(f,  X)  is,  in  general,  not  a 
single-valued  function  of  f.  A detailed  discussion  can  be  found  in  [3] 
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We  assume,  following  Barnhart  and  Street  [1],  that  the  particle 
number  density  n(f,  r)  consists  of  the  continental  and  maritime  compo- 
nents with  a variable  mix  ratio  so  that 


n(f,  r)  = unc(f,  r)  + (1  - u)  nM(f»  r)  (2) 

In  eq.  (2),  n^ff.  r)  and  nM(f,  r)  are  the  particle-number  densities  of 
the  continental  and  maritime  components  respectively.  The  quantity 
u(0<u<l)  is  the  mix  ratio. 


In  references  [1]  and  [2],  it  was  assumed  that  at  a given  rela- 
tive humidity,  the  rates  of  growth  of  all  the  aerosol  particles  were 
equal.  However,  because  of  their  different  chemical  compositions,  the 
growth  rate  of  a deliquescent  sea  salt  is  significantly  different  from 
that  of  the  condensation  nuclei  of  continental  origin.  In  fact,  it 
has  been  found  experimentally  [3]  that  y(f),  the  linear  mass-increase 
coefficient  of  an  aerosol  sample,  given  by 


U(f) 


1 - f 

m f 

o 


(3) 


where  mQ  is  the  mass  of  the  dry-state  aerosol  particles  and  n^,  the 
mass  of  water  uptake  of  the  sample,  can  be  very  different  for  different 
types  of  aerosols.  For  an  aerosol  particle  of  radius  rQ  in  the  dry 
state,  HSnel  [3]  has  obtained  the  following  relationship  giving  its 
radius  r at  relative  humidity  f: 


The  symbols  pq  and  p^  in  eq.  (4)  denote  the  density  of  the  dry-state 
aerosol  and  water  respectively.  The  second  term  on  the  right-hand  side 
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accounts  for  the  correction  due  to  the  finite  radius  of  the  particle. 

In  many  cases,  the  correction  term  can  be  neglected.  As  for  the  refrac- 
tive index  n(f.  A),  it  is  related  to  the  dry-state  refractive  index  of 
the  aerosol  n(f  = 0,  A)  and  n^(A)  of  water  by 


n(f.  A)  = nw(A)  + 


n(f  = 0,  A)  - nw(A) 


■](■ 


(5) 


3.0  CALCULATION  OF  THE  AEROSOL  EXTINCTION  COEFFICIENTS 


As  mentioned  in  the  previous  section,  the  aerosol  number  den- 
sity n(f,  r)  consists  of  two  components.  The  continental  and  maritime 
components  are  chosen  such  that  at  f = 0.7  they  are  the  same  as  those 
used  by  Barnhardt  and  Streete  in  [1] . By  means  of  numerical  fitting 
we  obtain 


nc(f  = 0.7,  r)  = 0.258  r'3'98 

(6) 

and 

nM(f  = 0.7,  r)  = 0.369  r"1  e‘°-748  r 

(7) 

Once  the  linear  mass-increase  coefficients  p(f)  are  given  for 
the  continental  and  maritime  components,  cqs.  (4),  (6)  and  (7)  enable 
us  to  calculate  n(f,  r)  for  other  values  of  the  relative  humidity  f. 

For  the  linear  mass-increase  coefficients  y(f),  we  have  chosen 
the  experimental  data  of  the  sea-spray  aerosol  [4]  i.e.  (Model  2 of  [3]) 
for  the  maritime  component  and  that  of  alpine  aerosol,  (Model  6 of  [3]) 
for  the  continental  components.  The  values  of  m(f)  we  used  are  in 
Table  I. 


UNCLASSIFIED 

5 


TABLE  I 


Linear  Mass-Increase  Coefficients 


f 

u(f) 

Maritime 

Continental 

0.5 

0.108 

0.032 

0.6 

0.136 

0.04 

0.7 

0.216 

0.0505 

0.8 

0.782 

0.14 

0.9 

0.617 

0.15 

0.95 

0.548 

0.105 

0.99 

0.468 

0.117 

At  f = 0.99,  the  uptake  of  water  is  46.33  times  the  mass  of 


I 


! 


the  dry-state  aerosol  for  the  maritime  component,  while  it  is  11.58  times 
for  the  continental  component.  Hence,  the  maritime  aerosol  particles 
grow  much  more  rapidly  than  the  continental  ones. 

Bamhardt  and  Streete  [1]  have  shown  that  to  calculate  8 in 
the  infrared  region  it  is  sufficient  to  consider  aerosol  particles 
with  radii  in  the  range  from  0.1  ym  to  20  ym.  In  this  size  range,  the 
correction  term  due  to  the  finite  radius  of  the  particles  in  eq.  (4)  is 
always  less  than  1%  for  both  the  maritime  and  continental  components. 
Thus,  this  term  will  be  neglected.  To  make  our  results  more  readily 
comparable  with  those  of  [1]  and  [2],  we  take  the  lower  limit  of  the 
integral  of  (1)  to  be  0.1  ym  for  all  values  of  f,  and  the  upper  limit 
equal  to  20  ym  when  measured  at  f = 0.7. 

Before  evaluating  the  efficiency  factor  for  extinction 

Qext  [a,  n(f , A)]  , it  is  necessary  to  have  the  complex  refractive 

index  n(f,  A).  The  quantities  nw,  n(f  = 0,  A)  entering  eq.  (5)  are 

summarized  in  Table  II.  In  this  table,  the  real  part  n ' of  the  index 

o 


UNCLASSIFIED 

6 


of  refraction  of  the  maritime  aerosol  component  in  the  dry  state  is 
taken  to  be  that  of  NACf  as  given  in  [5] . The  imaginary  part  is 
obtained  from  the  absorption  measurement  of  sea-salt  carried  out  by 
Volz  [6] . The  complex  refractive  index  for  the  continental  aerosol 
component  in  the  dry  state  is  that  used  in  [7] . The  index  of  refrac- 
tion of  water  is  taken  from  the  results  compiled  by  Hale  and  Querry 
[8].  Figures  la  to  Id  show  the  refractive  indices  of  the  continental 
and  maritime  aerosols  for  various  values  of  f as  obtained  from  eq.  (5), 
and  the  corresponding  dry-state  values . 


The  dry-state  densities  are  taken  to  be  pq (maritime)  = 
2.25  gm/c.c.  and  pq (continental)  = 1.8  gm/c.c. 
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FIGURE  1(a)  - Refractive  index  (real  part)  as  a function  of  wavelength 

for  maritime  aerosols.  The  number  accompanying  each  curve 
is  the  relative  humidity. 


FIGURE  1(b)  - Refractive  index  (imaginary  part)  as  a function  of  wavelength 
for  maritime  aerosols.  The  number  accompanying  each  curve 
is  the  relative  humidity. 
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FIGURE  1(c)  - Refractive  index  (real  part)  as  a function  of  wavelength 
for  continental  aerosols.  The  number  accompanying  each 
curve  is  the  relative  humidity. 


FIGURE  1(d)  - Refractive  index  (imaginary  part)  as  a function  of  wavelength 
for  continental  aerosols.  The  number  accompanying  each 
curve  is  the  relative  humidity. 
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Let  p„(f)  and  p..(f)  be  the  linear  mass-increase  coefficient 
L M 

of  the  continental  and  maritime  components  respectively  and  F^(f), 

FM(f)  be  the  corresponding  growth  factors  defined  by 

r = roFc(f)  (8) 

r - rQFM(f)  (9) 


Equation  (8)  describes  the  growth  of  an  aerosol  particle  of  the  continen- 
tal component  with  a dry-state  radius  rQ  and  eq.  (9)  is  the  correspond- 
ing relation  for  an  aerosol  particle  of  the  maritime  component.  Thus, 


PcCf) 


uc(f) 


CIO) 


and 

/ p o f \1/3 

FM(f)  = ( 1 + ^ ^ T J 

The  functions  Fc(f),  FM(f)  are  shown  in  Fig.  2 together  with  the  growth 
factor  used  in  [1].  From  the  assumption  that  no  aerosol  particles  are 
created  or  removed,  the  number  density  n(f,  r)  given  in  eq.  (2)  can  be 
related  to  nc(f  = 0.7,  r)  and  n^(f  = 0.7,  r)  as  follows: 


n(f.  r)  = u nc(f  = 0.7,  ~ r)  + (1  - u)  n^f  = 0.7, 

(12) 

The  calculations  of  the  efficiency  factors  of  extinction  can 
be  carried  out  using  the  Mie  theory  [9] . For  this  purpose,  we  have  used 
the  program  of  Dave  [10].  For  small  values  of  the  size  parameter  a, 
Rayleigh's  approximation  suffices.  We  have  used  this  approximation  for 
r < 0.5  pm.  The  integration  in  eq.  (1)  is  then  performed  by  means  of 
Simpson's  rule. 
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4.0  RESULTS  AND  DISCUSSIONS 

In  Figs. 3 and  4,  the  results  for  the  aerosol  extinction  coef- 
ficient 3 for  the  maritime  and  continental  components  are  shown  for 
various  values  of  relative  humidity  f.  For  the  maritime  component,  3 
is  a monotonic  decreasing  function  of  X in  the  wavelength  regions 
considered  for  f = 0.5  and  0.6.  When  f increases  beyond  0.7,  a broad 
minimum  occurs  around  11  ym. 

For  the  continental  component,  because  of  the  large  refractive 
index  of  the  dry-state  aerosol  around  9 urn,  the  extinction  coefficient 
shows  a marked  peak  in  this  spectral  region  at  low  relative  humidity. 

At  higher  relative  humidity,  this  feature  is  washed  out  with  the  accre- 
tion of  water  and  a broad  minimum  again  appears  around  10  urn. 

In  Figs.  5-11,  the  aerosol  extinction  coefficient  for  the  two- 
component  distribution  model  is  shown  for  various  mix  ratios  and  rela- 
tive humidities.  At  relative  humidities  equal  to  0.8  and  above,  3 
decreases  with  the  mix  ratio.  However,  for  relative  humidities  less 
than  0.8  this  trend  is  reversed  around  9 um.  This  is  due  to  the 
lower  index  of  refraction  of  the  continental  aerosols. 

Attempts  have  been  made  by  Barnhardt  and  Streete  [1]  and 
Hodges  [2]  to  compare  the  results  obtained  from  their  models  with  exper- 
imental values  of  the  aerosol  extinction  coefficients  measured  by  Yates 
and  Taylor  [11]  and  Kumick,  Zitter  and  Williams  [12]  . Unfortunately 
such  comparisons  cannot,  at  this  moment,  lead  to  definite  conclusions 
because  the  experimental  measurements  of  3 so  far  have  not  been 
accompanied  by  information  of  the  particle  number  distribution.  By 
changing  the  parameters  in  an  assumed  distribution,  one  can  obtain  a 
range  of  different  values  for  3.  In  fact,  Curcio  [13]  has  shown  that 
different  distributions  can  reproduce  the  same  extinction  coefficients. 
However,  it  is  of  interest  to  compare  the  present  results  with  those  in 
[1]  and  [2]  based  on  a simplified  model. 
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In  general,  the  values  of  3 obtained  in  [1]  are  higher.  This  is  because 
the  growth  factor  assumed  by  Barnhardt  and  Streete  is  larger  than  the 
empirical  values  we  have  adopted.  Take,  for  example,  the  case  where  the 
mix  ratio  is  0.5,  i.e.  the  continental  and  maritime  components  are  in 
1:1  ratio.  The  results  of  [1]  are  higher  than  those  obtained  here  for  all 
the  relative  humidities  except  the  highest  (f  = 0.99).  When  f = 0.99, 
the  agreement  is  reasonable  in  the  3-5  ym  region.  In  the  8-15  ym  region 
our  results  show  a dip  around  10  ym  which  is  absent  in  [1] . In  Fig.  2 
we  see  that  the  growth  factor  used  in  [1]  is  between  the  empirical  values 
for  the  continental  and  maritime  components  when  f = 0.99.  It  can  also 
be  noted  that  the  presence  of  a dip  in  the  values  of  8 in  the  10-12  ym 
spectral  region  is  characteristic  for  many  fogs  and  clouds  [14]  . 

In  [2] , Hodges  has  assumed  the  same  single  growth  factor  as 
in  [1]  for  both  the  continental  and  maritime  components.  He  has  also 
assumed  that  both  components  had  the  same  refractive  index  except  which 
now  could  be  complex.  The  imaginary  part  of  the  refractive  index  has 
been  taken  to  be  equal  to  that  of  water.  From  Table  II,  we  can  readily 
see  that  the  assumption  n"  = leads  to  an  overestimation  of  the  aerosol 
absorption  coefficient  which  becomes  more  serious  in  the  longer  wave- 
length region  (A  > 10  ym) . This  is  particularly  clear  for  the  maritime 
component.  In  fact,  the  values  of  8 obtained  in  [2]  for  the  3-5  ym 
region  are  substantially  the  same  as  those  in  the  8-15  ym  region.  From 
our  results,  8 is  smaller  in  the  latter  region. 
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FIGURE  3 - Extinction  coefficient  as  a function  of  wavelength  for  maritime 
aerosols.  The  number  accompanying  each  curve  is  the  relative 
humidity. 
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FIGURE  4 - Extinction  coefficient  as  a function  of  wavelength  for 
r continental  aerosols.  The  number  accompanying  each  curve  is 

the  relative  humidity. 
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FIGURE  5 - Extinction  coefficient  as  a function  of  wavelength  for 

different  mix  ratios.  The  number  accompanying  each  curve 
is  the  mix  ratio. 
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FIGURE  6 - Extinction  coefficient  as  a function  of  wavelength  for 

different  mix  ratios.  The  number  accompanying  each  curve 
is  the  mix  ratio. 
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FIGURE  7 - Extinction  coefficient  as  a function  of  wavelength  for 

different  mix  ratios.  The  number  accompanying  each  curve 
is  the  mix  ratio. 


FIGURE  8 - Extinction  coefficient  as  a function  of  wavelength  for 

different  mix  ratios.  The  number  accompanying  each  curve 
is  the  mix  ratio. 
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FIGURE  9 - Extinction  coefficient  as  a function  of  wavelength  for 

different  mix  ratios.  The  number  accompanying  each  curve 
is  the  mix  ratio. 
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FIGURE  10  - Extinction  coefficient  as  a function  of  wavelength  for 

different  mix  ratios.  The  number  accompanying  each  curve 
is  the  mix  ratio. 
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FIGURE  11  - Extinction  coefficient  as  a function  of  wavelength  for 

different  mix  ratios.  The  number  accompanying  each  curve 
is  the  mix  ratio. 
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5.0  CONCLUSIONS 

Based  on  empirical  results  of  the  linear  mass-increase  coef- 
ficients and  the  dry-state  aerosol  refractive  indices,  we  have  studied 
the  influence  of  increase  in  relative  humidity  on  aerosol  extinction. 

In  doing  so,  we  have  used  different  sets  of  optical  parameters  for 
the  maritime  and  continental  aerosols  as  required  by  the  difference  in 
their  chemical  compositions.  For  both  components,  the  values  of  8 
increase  by  a factor  of  more  than  10  throughout  the  3-5  and  8-15  pm 
window  regions  when  the  relative  humidity  is  increased  from  0.5  to  0.99. 
For  a two-component  distribution  model,  the  aerosol  extinction  coeffi- 
cient in  the  shorter  wavelength  window  is  generally  higher  than  that  in 
the  longer  wavelength  window.  For  a relative  humidity  higher  than  0.8, 
a dip  in  the  value  of  8 around  10-11  pm  is  also  noticed. 

In  contradistinction  to  the  results  previously  obtained  by 
using  simpler  models,  the  influence  of  relative  humidity  on  aerosol 
extinction  has  shown  to  be  more  dramatic  and  a better  aerosol  trans- 
mission condition  in  the  8-15  pm  region  is  also  indicated. 
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